Copy 2() of 195 ccpies 


IDA DOCUMENT D-1085 


VARIABILITIES IN THE NATURAL AND NUCLEAR 
ENDOATMOSPHERIC ENVIRONMENT 


Emest Bauer 


® : al} 4 IC 


a ee SoS 
MAY 18 1993 


a April 1992 
. 3 B 


see ; 
? * 


Prepared for 
Strategic Defense Initiative Organization 


° | 985 17 of with 


os a EEE 


INSTITUTE FOR DEFENSE ANALYSES 


1801 N. Beauregard Street. Alexandria. Virginia 22311-1772 


@ 
IDA Log No. HO 91-40677 


SEA tens 


IDA DOCUMENT D-1085 


VARIABILITIES IN THE NATURAL AND NUCLEAR 
ENDOATMOSPHERIC ENVIRONMENT 


Ermest Bauer 


April 1992 


Approved for public release: distribution unlimited. 


IDA 


INSTI fUTE FOR DEFENSE ANALYSES 


Contract MDA 903 89 C 0003 
Task T-R2-597.12 


Igk = 1/11 Ogk (1) 


where n = particle number density which ranges from 2.5 x 1025 particles/m3 at sea level to 
1 x 1019 particles/m3 at 100 km, and Og, ~ 6 x 10-!5 cm? is the gas-kinetic mean collision 
cross-section of an air molecule. Thus Ip, varies from 0.07 um at sea level to 20 cm at 
100 km. 


Re-entering missiles first interact with air molecules near 100 km, where the gas 
kinetic mean free path lp, is of the oreer of missile dimensions. The peak heating of a 
typical high-performance ICBM occurs between 15 and 25 km. 
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Figure 1. Phenomenology In the Lowest 100 km of the Earth's Atmosphere 


The weather (i.e., clouds and rain) is confined to the troposphere, roughly the 
lowest 10 km of the atmosphere. Most dust and smoke and a large fraction of the blast 
effects are confined to this region, which contains approximately 75 percent cf the mass of 
the earth's atmosphere. 


A nuclear fireball for a high-yield burst is formed below 70-80 km; Table 1 shows 
that bomb X-rays (which account for roughly 75 percent of the bomb's energy) are 


typically transmitted in the atmosphere only above these altitudes since there isn't enough 
atmospheric mass at higher altitudes. 


Tabie 1. Transmission of X-Rays, Neutrons, and Gamma Rays In the Atmosphere 


Transmission through 10-km horizontal path 
<<, 
Altitude Density 1 MeV 1 MeV 
(km) (kg/m3) 1 keV 3 keV 10 keV | neutrons gamma rays 


Extinction Coefficient is: (0.0.7 (0.1 
crr2/g) emg) 


In Figure 1 we also show how the bomb energy is carried.!. Above 80 km it is 
mainly transported by X-rays. Below 70 km most of the energy is carried as thermal 
energy of the fireball, while below ~ 20 km, where the air density is much greater than at 
higher altitudes, blast and shock become most important as a damage mechanism. Note 
that while neutrons and gamma rays account only for a small fraction of the bomb energy, 
yet they are important as a kill mechanism in the 50-70 km altitude range. 


1 For more detail, see Section 2.1. 


Referring back to Figure 1, there are lots of test data for heights of burst (HOB) 
below ~ 10 km, some data from 20-90 km, and very few data above 95 km, so that with 
increasing altitude the phenomenology depends increasingly on analysis rather than on test 
data. 


Table 1 shows the atmospheric transmission through a horizontal] (constant density) 
10-km path for 1, 3, and 10 keV X-rays, and for 1 MeV neutrons and gamma rays, at 
altitudes below 100 km. We see that X-rays are absorbed between 40 km (10 keV) and 
90 km (1 keV), while the neutrons and gamma rays are absorbed between 20 and 30 km. 


In an atmosphere at constant temperature T, the density p(z) falls off with 
increasing altitude z as 


P(Z) = (Zp) exp ~ (z — Z)/H (2) 


where the atmospheric scale height H is given by the expression H = kT/Mg. 
Representative values for H are 7 km at altitudes below 100 km where the temperature is 
~ 200-250 K; above 200 km altitude, where T ~ 700-1500 K and the atomic oxygen is 
dissociated so that the effective molecular weight is typically 18 (rather than 29), H lies in 
the range 30 - 70 km. 


In Section 2 we give a brief overview of the nuclear endo-atmospheric 
environment, pointing out the differences in phenomenology associated with the large 
difference in density (factor 105) between sea level and 100-km altitude. This is followed 
in Section 3 by a survey of variable aspects of the natural environment (mainly in the dense 
lower atmosphere), noting that clouds (which occur frequently) can totally obscure electro- 
optical sensors, as can rain, dust clouds, and smoke from fires. 
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2.0 NUCLEAR PHENOMENOLOGY BELOW 100 km 


2.1 INTRODUCTION 


There exist significant data on eight U.S. nuclear explosions between heights of 
burst of 20 and 400 km (two above 100 km, 6 between 20 and 95 kam), as opposed to more 
than 100 between the surface and 5-10 km. Thus the low-altitude data base is very much 
better than that at high altitudes, and the warnings of uncertainty in the high-altitude data 
base given in Bauer, 1990, are not so critical here. However: 

(a) As is pointed out in Figure 1 (above), the density falls off by a factor of 106 
between the surface and 100 km, and thus there are significant variations in 
phenomenology throughout the region in which there is a significant 
atmosphere. 

(b) Because all atmospheric nuclear tests were conducted in 1962 and earlier, the 
caveats about the lack of UV and LWIR data given in Bauer, 1990, still apply, 
but they tend not to be so serious because of the higher density of the 
atmosphere in which most of the UV and IR radiation are absorbed. Thus the 
atmosphere tends to radiate and absorb as a black body, so that the details of 
the radiating atomic and molecular species are normally not as critical as at 
higher altitudes where the spectra! variation in atomic and molecular radiation is 
critical. 

Much of the energy of a bomb is emitted initially as 1-10 keV X-radiation, which is 
absorbed in < 10 km of atmosphere at altitudes below 70-90 km (see Table 1). A fireball is 
produced at the relatively high densities corresponding to altitudes below 70-90 km by the 
energy emitted from the bomb--both X-ray energy (about 75 percent of the total energy) 
and bomb debris (20-25 percent of the total energy). The dimension of this fireball is of 
order 1 km for 1 Mt yield. At very low altitudes much of the fireball energy is dissipated 
by strong blast/shock effects, with the balance emitted as thermal radiation. If the burst is 
sufficiently close to the surface a large dust cloud is produced. 


A small fraction of the energy of a bomb is emitted as nuclear radiation. The 
neutrons and gamma rays--while they may account for only 0.1-1 percent of the total 
energy of the detonation--can have very significant effects (especially on sensors, 
electronics, etc.) because they penetrate the atmosphere very effectively. As is indicated in 


The hardness/vuinerability levels assumed are the following: 
* 10 cal/em? for X-rays 
¢ 1013 neutrons/em? 


* 50 cal/cm? for thermal (this is the incident fluence; if the surface reflectivity is 50 percent, 
then 25 cal/cm2 goes into the material at this range). 


* For blast/shock, we show the effective range for 1, 3, 10 psi, assuming the standard scaling, 
see, é.g., Glasstone and Dolan, 1977, p.100 ff. 
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Figure 2. Selected Effects Radil as a Function of Altitude for a 1 Mt Weapon 


* Above about 70 km, the damage effect of soft X-rays normally predominate. ! 


The term "fireball" is used--conventionally but inconsistently--to refer to two 
distinct concepts: 


* For altitudes below ~ 70 km, it is a visually defined volume in which 
essentially all the yield of the bomb is deposited. This volume of heated air is 
the source of the blast/shock wave and of the thermal radiation from the 
weapon. 


Above 100 km, the soft X-rays, which carry some 70-75 percent of the total bomb energy, are not 
absorbed in distances less than a few hundred kilometers. Thus exoatmospheric nuclear effects are 
produced mainly by the absorption of that 20-25 percent of the total yield that is carried by bomb debris 
as kinetic energy, which tends to be absorbed in the low-density air and is reradiated in the UV spectral 
range. 
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* The concept of a fireball is still useful for higher altitude bursts (for high 
yields, up to 150-200 km), but due to the lower ambient density and thus the 
longer mean free path for radiation, these higher altitude fireballs will have less 
well-defined edges, and contain 2 fraction of the total bomb yield that decreases 
with increasing alttude. 


2.2 PHENOMENOLOGY FOR BURSTS BETWEEN 20 AND 80 km: 
THE "HIGH-ENDOATMOSPHERIC REGIME" 


Reference to Figure 1 shows that the atmospheric density falls by a factor of 106 
between mean sea level and 100 km. Much of the energy of a bomb is emitted as X-rays: 
at altitudes below 70 km they are absorbed within 1-10 km or less, giving rise to a fireball. 
At the lower altitudes (below ~ 40 km for a 1 Mt burst) this fireball is smaller than the local 
atmospheric scale height? H = kT/Mg ~ 7 km and thus rises as a buoyant bubble, entraining 
outside air during its rise. At higher altitudes (above ~ 70 km for a 1 Mt burst) the ambient 
density is much lower and there is less entrainment, so that the fireball rises more rapidly, 
“ballistically,” and overshoots its final stabilization altitude.? Note that the upward speed 
of the fireball increases greatly with altitude. 


An air parcel of radius R is said to rise buoyantly if R << H, and ballistically if 
R >> H. Buoyant rise, which occurs at relatively low altitudes, corresponds to relatively 
slow, adiabatic rise of the air parcel, which maintains a uniform pressure; ballistic rise, 
which takes place at lower densities, is relatively rapid with pressure varying in the parcel. 
(See, e.g., Sowle, 1977, pp. 480 ff and 505 ff.) 


Table 3 indicates how the scale of phenomena changes with increasing altitude, and 
Figures 3 through 6 illustrate disturbed environments for three different modeling regimes 
as determined by burst altitude. Note both that the scale of the disturbed region increases 
as one goes up in altitude (and down in density) and the phenomenology changes. 


Regarding the phenomenology in the different altitude regimes, Figure 3 for a near- 
surface burst* comes from Glasstone, 1964, pp. 89-90. Note that while the fireball does 
not touch the ground (definition of an air burst), yet the afterwinds can sweep up a 
relatively small amount of dust.° 


2 Which is defined in Eq.(2), Section 1. 

3 At intermediate altitudes (40-70 km for 1 Mt) the behavior is intermediate between the “buoyant” and 
“ballistic” limits. 

This discussion of the physics is very useful, but some of the numbers are slightly inconsistent with 
current models. 


5 Cf. Figure 7b below; for both of these examples, the SHOB (“scaled height of burst”) is 
6,500 fyMt!/3, 
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Table 3. How the Scale of a Nuclear Detonation Changes with Altitude 


Altitude (km) 150 
Ambient density (kg/m?) ; : 1.8 x 10-9 


Gas-kinetic mean free path (m) 


Fireball size (km, for a 100 kt weapon) 
att= 1 sec 
at t = 30 sec 
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Figure 3. Chronological Development of a Near-Surface Burst 
(20 kt @ 0.5 km, 1 Mt @ 2 km) 
(Source: Glasstone, 1964) 
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Nuclear Explosion. (Source: Rausch, et al., 1988) 
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Time after Burst (hr) 


Fast cloud spreading (Curve IV of Fig. 10) 


Horizontal Cloud Width (km) 800 1000 Large extent 
Percent Degradation of Target Contrast 0 Rapid recovery 


Slow cloud-spreading (Curve II of Fig. 10) 


Horizontal Cloud Width (km) 130 Small extent 
Percent Degradation of Target Contrast >99 10 Stow recove 


b. impact of a 1-Mt Surface Burst 


Figure 11. Surveillance from Space Through a Nuclear Dust Cloud. 
(Source: Bauer, 1985) 
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Table 5. High and Total Cloudiness at Representative Locations 
in the Northern Hemisphere. 
@ (3DNEPH data from Malick and Allen, 1978, 1979) 


Location 


China Lake, CA 
Grand Forks, ND 
Maui, HI 

Hudson Bay 

N. Atlantic S 

N. Atlantic N 

Jan Mayen Is. 
Thule 

Barrow, AK 


Arabian Sea 
Teheran 

| lonian Sea 

| Moscow 

Tyuratam 

Lop Nor 

Viadivostok 
Japanese Trough 
Anadyr 

» Murmansk 


. Highfotal cloudiness means, for example, that at China Lake in January high clouds occur 
0.17 of the time and total cloudiness occurs 0,38 of the time. 

. These data come largely from downward viewing satellites such as NOAA-6 and DMSP, which 
tend to under-report optically thin clouds. 

. Clouds are reported as present when at least 1/10 of the appropriate field of view is covered 

» by clouds. 

. High clouds are those above 7 km, with the altitude determined by the effective radiative 
temperature in the 10- to 12-um infrared band as compared with the atmospheric temperature/ 
altitude profile. 

. High clouds are thus mainly moderately thick cirrus or cirrostratus, plus some cumulonimbus 
(thunderciouds) at the lower latitudes (<30°). 
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SOURCE WORGNICZ. 1972 
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Figure 12. Frequency of Seasonal and Annual Occurrence of Rain at Potsdam, 
Germany (Woronicz, 1972) and (Annual Only) In Washington, D.C. (Crane, 1980) 


In Fig. 12 we also show the annual mean precipitation figures for Washington, 
D.C., (from Crane, 1980) which differ somewhat from the Potsdam data by showing more 
high-intensity precipitation, as one would expect from the higher frequency of convective 
activity (thunderstorms, etc.) at lower latitudes. 


Crane, 1980, shows data comparable to Fig. 12 for a variety of climatological 
regions. When the rain rate exceeds 1 mm/hour, the zenith attenuation at 60-100 GHz 
exceeds 1.5-2.5 dB. 
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APPENDIX A 


REVISION OF DNA NUCLEAR CRATER 
SPECIFICATIONS 


Source: Frederickson, 1991 


APPENDIX A 
REVISION OF DNA NUCLEAR CRATER 
SPECIFICATIONS 


DNA has recently completed an "end-to-end" cratering validation program that 
resulted in dramatic reduction of the crater size thought to result from the surface detonation 
of modem strategic weapons. Although a major field exploration and several underground 
nuclear tests conducted in this program occupied the spotlight, numerical simulations were 
in many ways more central to DNA's success. This article recounts the integrated role of 
the numerical simulations, re-interpretation of existing nuclear data, and additional field 
events in the evolution of DNA's view on nuc! zar cratering. 


DNA developed a crater specification methodology for its 1972 Capability of 
Nuclear Weapons - Effects Manual Number 1 (EM-1) with the acknowledgment that the 
nuclear database was incomplete and probably inappropriate for application to strategic 
yield surface burst weapons. The cratering events conducted at the Nevada Test Site 
(NTS) employed low yield sources suspected to produce larger craters than modem 
weapons of strategic interest. Data from the several high yield cratering ev: :ts conducted 
at the Pacific Proving Grounds (PPG) were considered flawed by the atoll reef geology that 
was highly dissimilar to sites of interest. The 1972 EM-1 methodology was an attempt to 
reconcile these shortcomings. 


The strategic source surface burst crater specifications were based on high yield 
PPG data, calibrated to sites of interest by comparison of low yield nuclear and high 
explosive craters in various geologies. Figure A.1 depicts 1 Megaton crater profiles for 
two geology types as specified in 1972 EM-1. 
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Figure A.1. 1 Mt Contact Burst Crater Profiles for Two Generic Geologies 
as Specified by DNA EM-t (1972) 
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To address the database deficiencies, DNA developed computer code capabilities 
and applied them to numerically simulate cratering phenomenon. In the Benchmark 
Cratering Program (1976-1981), the 500 ton MIDDLE GUST high explosive event was 
conducted and its crater used to calibrate the numerical codes for simulation of a nuclear 
event on the same scaled geology. As seen in Figure A.2, the resulting simulated crater 
was markedly smaller and more bowl-shaped than the characteristically dish-shaped EM-1 
specifications. Simulations of high yield PPG events were also conducted and similar 


discrepancies with the reported crater profiles resulted. 
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Figure A.2. Comparison of 1 Mt Contact Profiles for Two Generic Geologies 
as Specifled by EM-1 (1972) and Predicted by DNA Benchmark 
Numerical Simulation (1981) 

Resolution of the discrepancy between crater specifications based on the existing 
but flawed database and the new numerical simulations became a central theme in DNA's 
Cratering and Ground Shock Program. The program sought to validate the simulation 
capability in separate, overlapping phenomenology components which, when placed end- 
to-end, spanned the entire nuclear cratering process. The phenomenology can be 
summarized in four Spica areas: 


eS arnnd aaediale Process occurs in first several SPATE for a 
Megaton yield event. 


* Conversion of coupled energy to ground motion field. This process is driven 
by the high pressure equation of state of ground materials and occurs in 
microsecond to several tens of milliseconds time regime. 


. n ion, transien r devel n ize. Thought to be 
dominated by the ejection of ground material from the ground, this period lasts 
100s of milliseconds for dry pcrous sites to several seconds for saturated soil 
Sites. 


“4 ; ff luti final f 
The key elements of the end-to-end validation process and the phenomenology 
component addressed by each element are identified in Figure A.3. 
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Figure A.3. Components of End-to-End Crater Validation Program 
Referenced to Phenomenology Regime Addressed 

The program addressed the last phenomenology components first. In the PEACE 
program (Pacific Enewetak Atoll Crater Exploration), DNA re-surveyed two high yield 
craters with the intent of determining whether small bow]-shaped ejecta/flow craters might 
have been formed within the overall reported crater dimensions. The theory to be tested 
was that such initial craters were subsequently altered by late time processes such as 
subsidence, slumping, or ocean washing that were not modeled in the simulations. The 
survey, conducted in 1983-1984, found compelling evidence that this was the case. For 
the first time, due to the insight gained from numerical predictions of the cratering process, 
a survey had been conducted that looked for the right data in the right places. Previous 
surveys had quantified crater extent based on observed deformations that had nothing to do 
with the environments of interest to vulnerability/survivability studies. With this new 
understanding of the Pacific craters, the large discrepancy with numerical simulations was 
gone. However, the unique atoll reef geology of the PPG meant this accomplishment was 
a necessary but not sufficient test for validation of simulations when applied to sites of 
strategic interest. 

Past tests provided craters and ground shock data in good agreement with pretest 
numerical simulations. The numerical simulations indicated that strategic yield sources 
would produce craters one-third to one-fifth the scaled size produced in these event: due to 
the relative inefficiency of the x-ray coupling process relative to hydrodynamic coupling. 

This early time x-ray coupling piece of the end-to-end validation was still missing. 
Numerical simulations were again utilized to determine what might be accomplished in a 


cavity using the higher yield source necessary to produce adequate x-ray output. It was 
found that a space-time window would exist in the same sized cavity used in MINI JADE 
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and MILLYARD (Figure A.4) such that the pertinent coupling physics could occur prior to 
the arrival of signals from the cavity walls. 


11 METER-RADIUS 
CAVITY 


Figure A.4. Contiguration of MIN! JADE and MILLYARD 
Underground Nuclear Cavity Tests 

Today, DNA relies on numerical cratering and ground shock simulations as key 
integral parts of its experimental program. They are the basis for cratering specifications 
for near-surface bursts in EM-1, 1991. Figure A.5 compares 1991 EM-1 craters on two 
geology types to the profiles perceived in 1972. This dramatic shift in perception is based 
on the compelling evidence obtained in the highly successful field program discussed in 
this article. The current DNA reliance on numerical simulations is a result of the 
recognition that they provided the motivation for this program, enabled the success of the 
field activities, and today provide the means to apply this text experience to specific 
strategic weapon and geology combinations of interest. 
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Figure A.5. Comparison of Current DNA Specification for 
1 Mt Crater Profiles with 1972 Specifications 
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